Non-Markovian continuous quantum measurement of retarded observables 
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We reconsider the non-Markovian time-continuous measurement of a Heisenberg observable x 
and show for the first time that it can be realized by an infinite set of entangled von Neumann 
detectors. The concept of continuous read-out is introduced and used to re-derive the non-Markovian 
stochastic Schrodinger equation. We can prove that, contrary to recent doubts, the resulting non- 
Markovian quantum trajectories are true single system trajectories and correspond to the continuous 
measurement of a retarded functional of 3:. However, the generic non-Markovian trajectories are 
mixed state trajectories. 

This version merges an Erratum [PRL, in print] with my Letter [PRL 100,080401 (2008)], some 
corrections follow directly from the criticism by Wiseman and Gambetta [2(|, further corrections 
restore the validity of my Letter. Contrary to my suggestion there, the give n continuous measure- 
ment schemes cannot yield pure state trajectories but mixed-state ones [2(j. Yet, it is possible to 
retain my claim that the NMSE (5) describes true time-continuous measurement - with delay and 
retrodiction. 



^— > ■ 



CO 

> 

On 
00 

m 



o 



X 



Time-continuous measurement in quantum mechanics 
had long been an open theoretical issue because of the 
peculiarity of single quantum measurement itself. The 
Markovian theory emerged twenty years ago [j], 0, 0] 
from foundational considerations. The requests in quan- 
tum optics (and elsewhere) triggered another, partly in- 
dependent, line of progress with expanding applications 
0] ■ So far the Markovian theory of continuous measure- 
ment has become completely understood while the gen- 
eral non-Markovian one has remained an open issue even 
conceptionally. 

Markovian time-continuous quantum measurement 
theory 0, Q includes the Markovian stochastic 
Schrodinger equation (MSSE) of the post-measurement 
state vector ipt, cf. 1], as correlated with the read-out xt 
of the detector system that measures a certain Heisenberg 
observable Xt . A formal extension for the non-Markovian 
(even relativistic) case was published in ref. [5]. This 
work calculated the asymptotic state V'oo [x] only, in 
function of the whole read-out {xt',t 6 (—00,00)}, and 
determined correctly its probability distribution func- 
tional Poo [a;]. It could not interpret intermediate con- 
ditional states because the concept of continuous read- 
out was missing. This incomplete non-Markovian con- 
tinuous measurement theory remained largely ignored, it 
has not been improved or advanced. Meanwhile, Strunz 
found non-Markovian quantum trajectories [f| and we in- 
vented their non-Markovian stochastic Schrodinger equa- 
tion (NMSSE) This NMSSE and its modifications 
have been studied in subsequent works 0, [l(J 11, [l2j ]. 
Like in the Markovian case, one expected that the solu- 
tions of the NMSSE turn out to be realizable on a sin- 
gle copy of our quantum system via infinite many von 
Neumann detectors coupled to it. Such realizability the- 
orem holds for the solutions (quantum trajectories) of 
all diffusive MSSE [Hj]. Yet, Gambetta and Wiseman 
conjectured that the solution of the NMSSE can not be 



observed on a single system [10( ; I wrote cautiously \U 
these non-Markovian trajectories_can not be realized by 
any known way of monitoring 



15] 



The present work reaches the positive conclusion: the 
non-Markovian trajectories are measurable single system 
trajectories. A particular example can be the contin- 
uous measurement of a Heisenberg coordinate Xt with 
detectors of finite inertial time 1/A. Then the measured 
quantity becomes, e.g.: 



% = A 



-A(t- 



T d<7 . 



(1) 



Our work includes the more general case, see eq. 
later. We describe the detector system and prove that 
the NMSSE is indeed the equation of the continuously 
measured state. The proofs are based on the approach of 
refs. 0,0,0,110,113, 

an independent direct proof might 
be subject of future research. The knowledge of the su- 
peroperator formalism is a request; it can be understood 
from [B|, [l6| or learned from 17 1 . 

Stochastic unraveling. Assume that a Heisenberg vari- 
able Xt of the system couples for times t > to a har- 
monic reservoir variable whose equilibrium correlation 
function a(r — a) will determine the reduced dynamics 
of the open system density operator p t : 



p t = M t po 



(2) 



For simplicity, let a(r — a) be real. Then the evolution 
superoperator M. t takes the following compact form [lj|: 

Alt=Tcxp^-i J c\ T J dax T:A a(T - a)x a ,A^ ■ 

(3) 

Superoperator notation x T:l \0 means [x Tl O] for any op- 
erator O standing to the right of x t ,a and T prescribes 
time-ordering for all Heisenberg (super) operators stand- 
ing to the right of T. 
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We could consider the reservoir as detector of Xt . Tech- 
nically, it is more tractable if we consider standard von 
Neumann detectors hence we replace the reservoir by 
them. However, we require that their influence on the 
system be the same as the reservoir's. We assume, for 
simplicity, that the detectors are able to fully monitor the 
system's trajectory ipt[x] for all time t > 0, in function 
of the detection read-out {x t ;t e [0,i]} whose probabil- 
ity distribution is denoted by p t [x] . Then the stochastic 
mean of the trajectories will reproduce the open system 
evolution: 



p t = MV> t [a#J[ 



(4) 



for all t > 0, since the detector's influence is the same as 
the reservoir's. We say that the trajectories 4>t[x] unravel 
the open system dynamics ©• 

In the Markovian special case a(r — a) = g 2 S(r — a). 
Then the conditional state vector ipt[x] satisfies the 
MSSE [J, B[3|. The NMSSE 0, i became a candidate 
of being the equation of non-Markovian continuous mea- 
surement of Xt ■ Here we use the simple real-noise version 
13, U, I2|- For the unnormalized state vector ^t[z], the 
NMSSE reads: 



d* t [z] 
dt 



z t x t ^t[z] - 2x t a(t 
Jo 



t)— dr , 



(5) 



where z r is a real random variable for r € [0, t]. The 
true post-measurement state is obtained via normaliza- 
tion ^t[z\ — ^t[z]/\^t[z\\\ . The probability distribution 
of z is the following: 



Pt[z]=G m [z] ||* t [z]|| 5 



(6) 



where G[ 0j t][z] is defined by (|30|) . With this statistics, the 
solutions tpt[z] unravel the non-Markovian open system 
dynamics (|2|3j) : 



pt = M1h[z]4[ 



(7) 



Although to calculate the analytic form ((6]) of p t [z] would 
be cumbersome, it follows from the method |8( that 



Mz t = 2 a(t- a)(x a ) t da 



(8) 



where {x a )t is x^s quantum expectation value at time 
t in the conditional state ipt[z]. This suggests that the 
NMSSE (O measures the retarded functional of xt rather 
than Xt itself. Compared to the Markovian case, there 
has been one serious issue left: Whether the trajectory 
ipt[z] can, like the Markovian trajectories, be realized on a 
single system by sequential von Neumann measurements 
of which z t is the read-out? We answer in the positive 
and construct the corresponding von Neumann detectors. 

Non-Markovian measurement device. The construc- 
tion will be very similar to the Markovian one 0, 18 1 in 



that we replace the reservoir by a dense sequence of stan- 
dard von Neumann detectors. To learn what happens, let 
us first consider a single von Neumann detector of initial 
density matrix Do(x;x') and couple it to our system at 
time r in order to measure the current Heisenberg op- 
erator x T . Following von Neumann (last three pages in 
[3]), we choose S(t — T)x T (—id/dx) for the interaction 
Hamiltonian. We can write the initial composite state 
of the detector+system as D (x;x')po- Fortunately, we 
can and shall restrict all forthcoming calculations on the 
elements x = x' since we shall eventually collapse on (or 
trace over) the pointer coordinates. After the interaction, 
the total state becomes entangled at r and the pointer x 
gets shifted by x T : 



D (x;x)p — ► D (x - x TiL ;x - x T ^)po 



(9) 



In superoperator notations x t ^O = x t O and x t ^rO = 
Ox T . It is the read-out of the pointer x that turns the to- 
tal state into the following conditional post-measurement 
state, depending on the read-out, of the system alone: 



p{x) 



1 



p(x) 



Dq(x - XrX] 



Xt,b)po 



(10) 



The read-out x has the probability distribution p(x) 
whose expression follows from the normalization of the 
above conditional state: 



p(x) = trDo(x — x T 



,l] x 



,R)Po 



(11) 



Now, let us choose a fine discretization r = ne of the 
time, n — 0, ±1, ±2, .... We install an infinite sequence 
of von Neumann detectors, they could be numbered by 
the integers n but we label them by the discretized times 
t = ne. The pointer coordinates of the detectors will 
be respectively denoted by x T . The detector of label 
t = ne measures the Heisenberg operator x T of the sys- 
tem via the mechanism (|9l [TTj) provided we switch the 
von Neumann interactions on. We do so for the non- 
negative labels, i.e., we choose the interaction Hamilto- 
nian J2 T >o ^ — T)x T {—id / dx T ) . 

We depart from the Markovian construction and as- 
sume initially correlated detectors. Let their initial wave 
function be: 



4>o [x] 



= VJ7 cxp ( — e 2 x t ol(t 



cr x a 



(12) 



where the summation extends for all discretized values 
of both t and a. The notation [x] anticipates the con- 
tinuous (or weak measurement) limit [1, 18| e — > where 
the above wave function becomes the square root of the 
Gaussian functional (f28|) . i.e.: <fio[x] — ^/G[x]. We carry 
out the explanation in the continuous limit. The total 
initial density matrix reads: 



Po[x;x'} = y/G[x]p y/G[x' 



(13) 
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As we switched on the detectors of labels r > only, 
at time t > each pointer coordinate x T with r S [0, t] 
will have been shifted by x T and the following composite 
state emerges [cf. ©]: 



p t [x\x] = T^Jg[x - 9 [0 ^x L ]^G[x - 9[ 0yt] x R }po , (14) 

where #[o,t] denotes the characteristic function 0[o,t](''") of 
the period [0,t]. This can be written into the following 
compact form: 



p t [x; x] = TG[x - 0[ 0>t ]X c }M t po 



(15) 



using the eqs. (|3|28|) and the superoperator notation 
x c O = O}. This remarkable novel form guarantees 
explicitly that the reduced density matrix pt of the sys- 
tem satisfies the open system evolution (|2l3p as it should. 
Indeed, the tracing over the detectors' Hilbert space is 
equivalent with the functional integration of the diago- 
nal elements (IT51) over all x T , which cancels the factor G 
and leaves us with ©. 

Continuous read-out. It is crucial to realize that the 
true time-evolution of the system's conditional state de- 
pends on our chosen schedule of reading out the pointers 
x T . We can read out any x T at any time since all detec- 
tors are always available. Of course, we better read out 
the value x T at a time which is later than the label t of 
the detector because the detector will only have coupled 
to the system at time r. Hence, a natural schedule is that 
we read out x T immediately at time r. Hence, until any 
given time t > we would read out all pointers x T for the 
period [0, t] and no others. To calculate the conditional 
post-measurement state pt[%] of the system at time i, we 
trace (integrate) the total density matrix (fl5|) over all x T 
with t <£ [0,t]: 



1 



Pt[x\ 



Pt[x\ 



p t [x;x\ 



IT dXr 

rj£[0,t] 



(16) 



This post-measurement density matrix pt[x] of the sys- 
tem depends on the read-outs x T of t from [0, t] only. By 
substituting (fl5|) . we obtain: 



1 



Pt[x\ 



Pt[x\ 



}M t po 



(17) 



where G[ t ][a;] is the marginal distribution of G[x], simi- 
larly to (|30p . This is our ultimate equation for the non- 
Markovian continuous measurement of the Heisenberg 
observable Xt , completing the theory [f| (which only gave 
Poofx]). Recall that, as always, the denominator p t [x] as- 
sures trp t [x] = 1 as well as it yields the probability dis- 
tribution of the read-outs. Contrary to our assumption, 
the state (17) is not pure even if it started from a pure 
Po; the continuous readout of Xt cannot provide sufficient 
information for a pure state ipt [x] , as shown by Wiseman 
and Gambetta [20(. 



In order to find the measurement process that corre- 
sponds to the NMSSE ([5]), we alter our read-out sched- 
ule. Instead of the Heisenberg variables {x t ;t £ [0,t]} 
we read out the following linear functional of them: 



/oo 
a(r — a)x a da 
-oo 



(18) 



which we also write as z = 2ax. We re-express the total 
density matrix (|15|) in the new pointer variables: 

Pt[z;z] = TG[z - 2a9 [0yt] x c }MtPa , (19) 

where we used the identity G[x] = Jacobian x G[z]. 
Again, we suppose that we read out each pointer of label 
t (i.e.: z T ) at time r. Until time t > 0, this schedule 
implies that all pointers z T for the period [0, t] are read 
out and the rest of them are not. The conditional state 
of the system is defined by: 



Pt [z\ 



T£[0,t] 



(20) 



which transforms (JTUJ) into 

~ r , 1 



Pt[z\ 



Pt[z_ 



T G m [z - 2a9 [0A x c ]M t p , (21) 



where Gr 0it i[z] is the marginal distribution (|30[) of G[z\. 
This is our ultimate equation for the non-Markovian con- 
tinuous measurement of the observable 



Zt 



2 / a(t — o-)x a da 
Jo 



(22) 



which is a retarded functional of the Heisenberg variable 
x r . This interpretation of pt [z] can shortly be inspected. 
Recall that at time t we read out the pointer of label t, 
i.e.: Zt- The factor Gro,t][^ — 2a9[ Q ^x c ] in the expres- 
sion (12111 of the measured state shows that at time t the 
pointer z t localizes around (i.e.: measures) the observ- 
able (j2"2")l . The eq. © holds between the read-out z t in 
pTjl and the retarded variable z t (|22[) ; instead of the di- 
rect proof we are going to prove the complete equivalence 
of the NMSSE © with our construction summarized by 
eq. (ED)- 

Stochastic Schrodinger Equation. We are going to 
prove that the NMSSE ([5]) governs the evolution (j2Tj) . 
Let us find pt[z] in the form: 



^ = ■zj- ] G m[z]^t[z}^t[z] 

Pt[ z \ 



(23) 



where ^t[z] is the unnormalized conditional state vector 
of the system. Taking the trace of both sides, the norm 
condition yields exactly the p t [z] ([6]) that belongs to the 
NMSSE ©. Inserting $F5§ as well as p = tp tpl into 
(|2"Tj) . it reduces to: 

*t^]* t f [z] = -, 1 r T G m [z - 2a9 [0it] x c }M t 4>o4>l ■ 



G[o,t] [z] 



(24) 
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Substituting eqs. ([3]) and (|30|) . the r.h.s. factorizes and 
we can write equivalently: 



\&t [z] — Texp I / z T x T dr — / dr / dax T 
\Jo Jo Jo 



a(r - cr)2 o .)V>0 



(25) 

This is the solution of the NMSSE as can be 

seen by substitution. That completes our proof. 

Delayed continuous readout. Unfortunately, the chosen 
readout schedule alters the reduced dynamics (2) because 
the detector modes (18) are not retarded, hence the cou- 
pling between the system and the detector mode z T con- 
tinues after time r, cf. Ref. [jo| . It ceases, nonetheless, 
at r + T provided T > is much larger than the reser- 
voir correlation time so that a(T) = be already a good 
approximation. We can thus keep the reduced dynamics 
(2) invariant if we apply continuous readout with a finite 
delay T . We read out each pointer of label r (i.e.: z T ) at 
time t + T . The conditional state (20) of the system at 
time t > T must be replaced by: 



p t [z; delay^T] 



1 



Pt[z;z] Y[ dz T 

T$[Q,t-T\ 



pt[z; delayer] 

(26) 

It turns out that pt[z; delay=T] = p t ~T[z], i.e., the statis- 
tics of delayed continuous readouts obtained until time 
t is identical to the statistics of zero-delay (and all-in- 
one [2(|) readouts until time t — T. The delayed-readout 
state obviously coincides with the following average of 
the zero-delay- readout states (20): 



p t [z; delay =T] 



1 



Pt-r[z] 
1 

Pt-r[z] 



Pt[z] Pt [z]\{dz T (27) 

T£[t-T,t] 

1p t [z}lpt[z] p t [z] Y[ dz T . 

re[t-T,t] 



The second equality follows from the insertion of p t [z] = 
ifj t [z]tpl[z] where, according to Eq. (23), the pure 
state ipt[z] must be the normalized solution ipt[z] = 
*t[z]/||* t [z]|| of the NMSE (5). As we see, ip t [z] does 
not directly describe a continuously measured quantum 
trajectory because the values z T for re [t — T, t] would 
belong to the all-in-one measurement at time t. Still, 
the above partial average of the pure states 4>t[z) over 
those fictitious z T does fully describe our (delayed) non- 
Mar kovian continuous measurement. The normalized so- 
lutions i\>t [z] of the NMSE (5) do correspond to retrod- 
icted pure states of the system, the proof and physical 
interpretation will be given elsewhere. 

Summary. We proved for the first time that both the 
formalism |5j of non-Markovian measurement theory and 
the NMSSE [?} are equivalent with using correlated von 
Neumann detectors in the weak-measurement continu- 
ous limit, i.e., with the continuous read-out of the val- 
ues of a given retarded functional of a Heisenberg vari- 



able on a single quantum system. Our merit is the con- 
structive proof of existence of the underlying standard 
quantum mechanical measurement process. The results 
should be generalized in various directions. We can in- 
terpret complex reservoir correlation functions, too, if wc 
include the mechanism of feed-back [3 . We might retain 
the original reservoir as detector [10j, to extract infor- 
mation by measuring the reservoir but without altering 
the non-Markovian reduced dynamics of the monitored 
system. Then the measured retarded observable might 
be identified by a reservoir field. (Theories advocating 
non-Markovian stochastic modification of quantum the- 
ory II, 12, [1H refuse the measurement interpretation of 



the stochastic field.) The concept of relativistically in- 
variant continuous measurement Q can be reconsidered 
for the intermediate states ifit [x] as well. Our work might 
lead to efficient numeric simulation algorithms or, con- 
versely, might make us understand why they don't exist. 

Appendix.- Let x T be a random time-dependent real 
variable and consider the normalized Gaussian distribu- 
tion functional of {x T ; r G (— oo, oo)}: 



G[x] = TVexp I —2 / dr/ dax T a(r — a)x a I 

\ J — oo J — oo / 



(28) 



a(r — a) is a real positive definite kernel, We define its 
inverse by a _1 (r— s)a(s— a)ds — 5(t— a). Introduce 
the normalized functional Fourier transform of G[x], too: 



G[z] 



~ ( 1 f°° f°° 
Af exp I — — / dr/ doz T a~ 1 ( 

\ Z' J — oq J —oo 



a)z a 



We need certain marginal distributions as well, e.g. 



G m [z]= f G[z] 11 dz T , 
rj£[0,t] 



(29) 



(30) 



and similarly for G[ 0; t][a:]. These marginal distributions 
are also Gaussian, e.g.: 



G[o, t ][z] =N[o,t] ex P(^2 J dT J daz T a [ ]t]( T ' (7 ) z <* 

. (31) 

where the restricted new kernel a^^r, a) is defined by 

Jo a [c) 1 t]( T ' s ) a ( s ~ <T )ds = S(t — a) for all r, a € [0, t\. 
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